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Abstract 

We study the rotational and translational kinetics of massive particulates (dust grains) ab- 
sorbing the ambient gas. Equations for microscopic phase densities are deduced resulting in the 
Fokker-Planck equation for the dust component. It is shown that although there is no stationary 
distribution, the translational and rotational temperatures of dust tend to certain values, which 
differ from the temperature of the ambient gas. The influence of the inner structure of grains on 
rotational kinetics is also discussed. 
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I. INTRODUCTION 



Recently there was growing interest to the investigation of composite media, called dusty 
plasmas, consisting of aerosol particles in a gas discharge. Besides numerous industrial 
applications dusty plasmas provide ample opportunities to studying in situ phase transitions 
in the dust component, intergrain interactions, grain charging etc. 

One of the peculiar features of dusty plasmas is that the average kinetic energy of the 
dust component, i.e., its translational temperature, may be considerably higher than the 
temperature of the ambient plasma. In details the problem was studied in the recent papers 
by Quinn and Goree |I|, 0, where in parallel with the experiments a model of the Brownian 
motion explaining anomalous kinetic temperature was developed. 

The kinetic description of dusty plasma was discussed in numerous theoretical studies 
||, f|, [|, §, [7|, ||, |§ . Generally, there are two ways plasma particles interact with dust grains: 
first, the scattering of a particle by grain electric field and, second, the direct impact of a 
particle on a grain surface. The latter process results in grain charging due to the higher 
mobility of electrons, it may change grain mass, heat its surface etc. In other words, as it is 
well-understood nowadays, the adequate statistical description of the dust component should 
take into account inner degrees of freedom, the most important among which is the grain 
charge. Kinetic consideration of charging process shows that absorption of small plasma 
particles by grains can result in inequality of the grain temperature and the temperatures 
of the light components even for the case of equal temperatures of electrons and ions[[|. |4[. 
To avoid confusion it should be noted that since the system is open there is no conflict of 
the latter result with thermodynamics. 

Under the conditions of experiments on fine grain synthesis |L0|, [□]] or etching [12|] grain 
mass should also be treated as a dynamic variable |13[ |14|, |15| . Thermophoretic and 



radiometric |L7j forces provided by inhomogeneous heating of the grain surface may also 
play an important role, under microgravity especially fTSf , that may necessitate inclusion of 
the temperature distribution inside grains into the kinetic description. Recent experiments 



on rod-like particulates |]T9| demonstrated complicated rotational motion of dust grains. In 
more general context, the rotation is also of importance for Brownian motion of the particles 
with energy supply [^(J. Although the rotation of spherical grains is hardly observable 
experimentally, it is of interest since it may influence translational motion and heating of 
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the grain surface. 

Being not quite complete, the above list of various process indicates that the dimension 
of the phase space required for the kinetic description of the dust component may be rather 
large: it tends to a rather large value of the order of twenty. 

The main purpose of the present paper is to develop the kinetic theory of the aerosol 
component taking into account the mass growth and the rotation of grains. Although our 
main impetus is dusty plasma, here we ignore the process of grain charging. This may 
be justified assuming the size of grains exceeds the Debye length of the ambient plasma. 
In other words, we treat the ion component as a neutral gas and neglect the influence of 
electrons. Another reason for accepting this simplification is that we are able to study 
the problem both analytically and by molecular dynamics simulation: computations with 
charged moving grains are on the brink of or beyond contemporary computer capabilities. 

Thus, we adopt here a following toy model of the aerosol component absorbing the am- 
bient gas. The dust component consists of spherical rotating grains with variable mass and, 
consequently, size and moment of inertia. Every atom hitting the grain surface is absorbed 
by it, transferring therefore, its momentum, changing the mass of the grain and its angular 
velocity. The process is inelastic since a part of projectile atom energy is spend for heating 
the grain surface. It is assumed that the size of grains is small compared to the mean free 
path of the ambient gas, however the gas distribution generally depends on the dust com- 
ponent. Our main finding is that although there is no stationary state of this system, the 
average kinetic and rotational energies of dust eventually tend to certain fixed values, which 
differ from each other and the temperature of the ambient gas. 

The paper is organized as follows. In Sec. II A we discuss the process of the elementary 
collision. Microscopic phase densities and corresponding equations generalizing Newtonian 
dynamics are introduced in Sec. |IIB| , the latter are used in Sec. |II C| to derive the Fokker- 
Planck equation (|20|) describing the evolution of the dust distribution. In Sec. |IID| we 
obtain the homogeneous solution of the Fokker-Planck equation and evaluate the effective 
temperatures. The numeric algorithms are described in Sec. [Ill A[ then we discuss various 
laws of collision used for simulations (Sec. [Ti FB]). The results of simulations, which are in 
fairly good agreement with analytical theory, are summarized in Sec. |111 Q 
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II. ANALYTIC THEORY 



A. Elementary collision 

We consider the dust component consisting of a number of spherical grains of variable 
masses, M, and of mass dependent radii, a(M). Since we are interested in both translational 
and rotational degrees of freedom, the state of a grain is described by the ten-dimensional 
vector, T = (R, P, G, M), where P is the grain's linear momentum and G is the angular 
momentum relative its center of inertia, R. The angles describing the rotation are irrelevant 
due to the sphericity of grains. The only collision process taken into account is the absorption 
of an ambient gas by grains, i.e., every atom colliding a grain is assumed to attach to the 
grain surface transferring its linear momentum, angular momentum and mass. 

We assume that the process of collision elapses in two stages. At the first stage the 
atom attaches to the grain surface. Let r and p be the projectile atom coordinate and, 
respectively, the momentum at the instant of collision; evidently, \r — R\ = a(M). The net 
angular momentum of the projectile atom and the grain prior to the collision is written as 
M = rxp + RxP + G, while after the collision M = R x P' + G'. Since after the 
attachment the center of inertia of the grain changes 

= mr + MR 

M + m v ' 

and P' = P + p, the conservation of angular momentum requires that 

Q , = G+ (r-R)x (Mp - mP) 

M + m 

It should be noted that both G and G' are independent of the reference frame. 

In what follows we combine the above expressions in the convenient short-hand notation 
for the process of collision: 

r^r- = A(r,r,p) = (=^«, (2) 

M+m 

„ (r-R)x (Mp - mP) „ , 

P + P,G+ K - ' \ y J -,M + m). 

M+m 

The leap of the center of inertia (|l|) results in the non-conservative character of the 
mapping ([|): one can easily verify that its Jacobian 
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det 



dA(T) 



M 3 



(3) 



dr 



(M + m) 3 



is less than unity. 

Strictly speaking, the attachment of the atom to the grain surface makes the new com- 
posite grain non-spherical. To avoid complications with non-spherical grains we assume that 
at the second stage of the collision some inner forces redistribute the grain mass in such a 
way that the grain shapes into a sphere. Evidently, this does not alter the grain phase 
variable, T. Although this assumption seems a fairly natural simplification, it may result in 
non-physical behaviour because it leaves out the energy required for the mass redistribution. 

Evidently, the process (|2]) is inelastic, that is, a part of the net mechanical energy is spent 
for heating the grain surface, which evidently should be positive. (The grain heating is of 
importance if one is going to consider the aerosol processes like thermophoresis.) Abandoning 
for a moment the presumption of sphericity, the amount of heat released in the asymmetric 
grain is written as 



where v = p/m — P/M is the atom velocity in the grain reference frame, m = mM/ (M + m) 
is the reduced mass, p = r — R, Iij and Zy = Z^ + fh (p 2 (% — PiPj) are the grain tensors 
of inertia prior and, respectively, after the collision. Expanding @ in powers of the small 
mass ratio, m/M, we obtain that AE — (v — f2 X p) 2 /2m > 0, where = Zy Gj is the 
grain angular velocity. In other words, the energy conservation does not debar the mass 
absorption. 

On the other hand, assuming the grain shapes into a sphere, the tensor of inertia is 
isotropic and depends on the grain mass only, Zy = Z(M)<5y, Zy = Z(M + m)5ij, and 
analyzing Eq. (|4]) one can find that AE is always positive if dI(M)/dM > a 2 (M). However, 
it is easy to verify for an arbitrary spherically-symmetric mass distribution that the latter 
inequality never holds. This means that some energy is required for the grain to shape 
into a sphere and there exist a certain range of impact parameters, such that the available 



AE = - {mv 2 + ClM;, 





5 



mechanical energy is insufficient for the mass redistribution. It should be noted that similar 
phenomena are also known in nuclear physics (e.g., ||21|| ). 

Simple estimation shows that the assumption of sphericity results in the violation of the 
energy conservation for atoms moving with characteristic velocities of the order of afl. In 
the following we assume that the energy of rotation is of the order of the gas temperature 
and the mass ratio, m/M, is small, i.e. the discussed effect is relevant for a very small 
group of projectile atoms with velocity less than ym/ 'Mvt, where vt is the gas thermal 
velocity. We ignore the influence of deviation from sphericity for this reason. Nonetheless, 
the discussed effect may play an important role for small clusters and, perhaps, it explains 
the complex structure of fine grains grown in a gas discharge. 

B. Microscopic phase densities 

In order to obtain the desired kinetic equation describing the grain motion we generalize 
the well-known Klimontovich approach ||22|| . To describe the whole system the microscopic 
phase densities are introduced: 



N d (r,t) = £)*(r-r a (f)), (5) 

a 

N n (p,r,t) = ^5(p-p^))5(r-n(i)), (6) 



where the indices a and % are used to enumerate grains and atoms, respectively. 
The equations governing the microscopic phase densities are written as 

dN d (T,t) fd P d „ d 
= \ — + —— + F d — 



dt \dt MdR dP 



+ K d -^jN d (r,t),= I d (T,R) (7) 



dNJp.r.t) ( d p d „ d , A . 

= ( 1 - + - ir + F n —)N n (p,r,t) 



dt \dt mdr dp / 

= In(p,r). ' (8) 

Here F d and F n are external forces acting upon grains and neutral atoms and K d is the 
external torque. The collision terms in the right-hand sides of Eqs. account for gas 
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absorption by dust grains; similar terms were introduced in |4j]. The convenient short-hand 
forms of these terms are 



I d (T, R) = J dp dr dT'N n (p, r)N d (V, R)a ( r -R,^-^, M' 

x[6(T-A(T',r,p))-6(T-T')], (9) 
UP, r) = - J dY a (r - R, £ - ^, m) N n (p, r)N d (T), (10) 

where the function A(T,r,p) is given by Eq. (0). The effective cross-section a(r,v,M) 
introduced here is 



a(r, v, M) = -2(r • v)6(-r • v)S{a 2 {M) - r 2 ) (11) 

and 8(x) is the Heaviside's step function. 

The collision terms (|§],[n|) are arranged in a following way. Suppose there is a trajectory 
described by a time-dependent radius vector r(t). If the trajectory crosses the sphere of 
radius a at the instant to, i.e., r(to) = a{M), then 

a(r,r,M) = 5(t-ta). (12) 

The step function in ( |TTj ) guarantees that the incoming intersection point, such that r(to) ■ 
r(to) < 0, is only taken into account. 

Substituting the definitions of the microscopic phase densities (fUD to Eqs. (0,0) and 
making use of ( |12"D one can verify that the collision terms are proportional to the sum of 
S (t — Tia), where Ti a is the instant of collision of the z-th atom with the a-th grain. At the 
time intervals between the collisions Eqs. (|7||f) describe plain Newtonian dynamics. However, 
at the instant of collision the microscopic phase densities change abruptly: 

N n (p, r, r ia + 0) - N n (p, r, r ia - 0) = -5(p - Pi)5(r - n), (13) 

that is, the z-th atom annihilates. Simultaneously, the a-th grain changes its position in the 
phase space as prescribed by the conservation laws ([|) 



N d (T,r ia + 0) -N d (T,r ia -0) = 

5(r-A(r a ,n, Pi ))-5(r-r a ). (14) 
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C. Fokker-Planck equation 

The purpose of this section is to expand Eqs. (pj,fL0|) in powers of the small mass ratio, 
m/M. Averaging over the ensemble and ignoring binary correlations, we may treat N d (T, t) 
and N n (p,r,t) as smooth one-particle distribution functions; in details the procedure was 
discussed in ||. Then we integrate over V in Eq. (jfty and take into account Eq. ([|) that 
results in 



I d (T) = [ dpdr I a ( r, — — — v, M — m) N d (r — —r, P — p,G — mr x v, M — m) 
J { \ M — m J \ M / 

{ M — m \ 1 

xN n [p,R+ —^r - a(r, v, M)N d (T)N n (p } R + r)\, (15) 



where 

m M v 1 

is the relative velocity. Evaluating Eq. ([T5|) we have also performed the change of variables 
in the integrand (Jffy: r — > R+ M ^™ r in the first term and r — > R + r in the second term. 

The next step is the expansion of Eq. ([T^) in powers of the grain radius, a(M), and mass 
ratio, e = m/M. Assuming that p/P ~ mav/G ~ e 1 / 2 and making use of the integrals 



J dra(r,v,M) = na 2 (M)v, (17) 

J dra(r,v,M)r = -^ixa\M)v, (18) 

J dra(r, v, M) ri r 3 = ^a A (M)v (<% + ^) , (19) 
we finally arrive at the Fokker-Planck equation for N d (T) 



dN d (T) d f dN d (T) \ d ( dN d (T) 

—*r = m \- s * Nd{T) + K -^p-\ + do; {n-dG- 

-fiN d {T) - e ljk a k —- U tJ } — — , (20) 



dRi {' v ' J dGj J dPj J dM 
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where is the unit antisymmetric tensor. Here the following kinetic coefficients are intro- 
duced: 



J(T) 
*(T) 

KijiT) 

7,(r) 
*(r) 

ny(T) 



vra 2 (M)m y dpvN n (p,R), 

na 2 (M) j dp vpiN n (p, R), 

^7ra 2 (M) | dpvp lPj N n (p,R) 
2m 



3M 

2 

m 



m 
4M 
2m 
3M 



vra 3 (M) y dpviN n {p, R), 
ra 4 (M) y dp v (^t; 2 - u^) N n (p, R), 

ixa\M) J dpv Vl N n (p,R), 
7ra 3 (M) / dpv lPj N n {p,R) 



(21) 



and t> is given by Eq. flU). 

The physical meaning of the most of the coefficients (|2lj ) and corresponding terms in 



Eq. (20) is fairly obvious. J is the mass flow at the grain surface, the last term in Eq. ( |20| ) 
provides the mass growth of the dust component. The coefficient is the drag force acting 
upon a grain, the quantities Kij and 77,-,- characterizes the diffusion in the momentum space. 
The term proportional to II^ is just the Archimedean force in a non-uniform gas. 

Since in the process of collision the angular momentum transferred to the grain is inde- 
pendent of its angular velocity, there is no drag torque analogous to the first term in Eq. (pop. 
This is the evident consequence of the adopted model. The drag torque may appear if one 
takes into account the non-sphericity of grains or inelastic scattering of atoms by the grain 
surface. Within the present model the drag torque is provided by spatial gradients and it 
is characterized by the coefficient ov Since in the low velocity limit, P/M <ti p/m, Oi is 
proportional to s», the drag torque may appear if the curl of the drag force is non-zero. 

Of interest is the term proportional to 7^ in Eq. (pup. It is non- vanishing if the grain 
moves relative to the ambient gas and arises due to the change of the center of inertia in the 
process of collision (fj) or, in other words, 7« describes the migration of the center of inertia 
due to the asymmetric bombardment of the grain by gas atoms. Comparing the left-hand 
side of (0) with this term we see that % appears as an addition to the grain velocity, P/M. It 
often happens that the ambient medium impose some forces at the particulate, an example 
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is the drag force, S{. In this case one may say that the medium modifies the second Newton's 
law. Here we face with the example of the first Newton's law altered by the ambient gas. 

By the order of magnitude the drag force, Sj, is proportional to e 1//2 , the effective torque, 
a", ~ e 3 / 2 , all other coefficients are of the order of e. 

In the lowest-order approximation the kinetics of the gas component is reduced to the 
absorption. The corresponding kinetic equation is readily obtained from Eq. (|TT}): 

=~ fdT rra 2 (M)^N n (p)N d (T) + I a (p, r), (22) 
at J m 

where I a (p,r) stands for other dissipative processes, which were excluded from the above 

derivation. These may be, e.g., interatomic collisions, gas creation by external source etc. 



D. Effective temperature 

Suppose that the gas distribution, N n (p), is homogeneous and isotropic. Since we as- 
sume that the grain velocity is small compared to the gas thermal velocity, we neglect the 
momentum dependence in all kinetic coefficients ([H]) but the drag force, Sj, for which the 
first order term of expansion in powers of P/M should be kept. Therefore, there are only 
three non-zero kinetic coefficients which are expressed in terms of mass flow, J, and 
normalized energy, a: 



J 



P 



where 



6mJ 

The Fokker-Planck equation fl20D is now reduced to 



3M 

Kij = dijJa (23) 
Vij = 5 ij ^a 2 (M)Ja, 



« = ^^ / dpp 3 N n (p). (21) 



dN d (T) 
dt 



d 

m 



J dN d (T) 
PiN d (F) + aJ 



+ -aa\M)J 



3M y ' dPi 
d 2 N d {T) dJN d {T) 



dddGi 



dM 



(25) 
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Suppose there is no dispersion over the grain mass. Then we seek the solution to Eq. fl25|) 
in the form of 



N d (T) = S(M-ii(t))f(P,G,ii(t)). 



(26) 



Evidently, if all the grains are of the same mass and of the same mass growth rate, we 
are able to use the current value of mass, fi(t) for parameterizing the temporal evolution in 
the phase space. Substituting the latter ansatz to Eq. (|25| ) we find that 



dt 
<9/i 



J 



— ^M- + <*A P / + -« 2 M«A G /, 
3/i dPi J 2 t " ,i 



(27) 
(28) 



where A p and Ac stand for the Laplacian operators acting on the corresponding variables. 

It should be noted here that the normalized energy, a, is formally independent of mass, 
M. However, it may depend on /i due to the possible time variation of the atom distribution, 
N n (p,t). 

Eq. (p8|) is reduced to the diffusion equation by changing the variable P — > x = fj}' 3 P. 
The latter is readily solved resulting in 

MP'.G') 



f(P,G,n) = / dP'dG' 



x exp 



(2vr)3(A; 1 (/i)fc 2 ( /U ))3/2 
{^P-^P'f {G-G'f 



2fci(/x) 



(29) 



where 



a(yu)/x 2 ^ 3 (i/i, 



mo 



a(fi)a 2 (fi)dfi, 



fio is the value of mass at t — and fo{P, G) is the corresponding initial distribution. 

As it follows from Eq. (^9l), f(P, G, fi) eventually tends to the Maxwellian distribution. In 
order to obtain the parameters of asymptotic distribution one have to evaluate the momenta 
of f(P,G,fi). First, it is readily checked that the dust density, rid = f dPdGf(P,G, fx), 
is independent of fi. Then, we evaluate the linear momentum dispersion, A t (fi), and the 
angular momentum dispersion, A r (/i): 
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3n d J 



I A *(^o) + ( 3 °) 

A,(/i) = ^- / dPdGG 2 f(P,G,fx) 

= A r (/i Q ) + A; 2 (//). (31) 



Explicitly, the asymptotic distribution is given by 



foo{P,G,fJL) 



(2^)3(^(^)^(^)3/2 



"^K"^)"^))- (32) 

It is a matter of direct substitution to verify that the latter distribution satisfy Eq. (|25|) 



subject to relations (|30| , |3T 



Now we introduce the translational temperature, T t (fi) = A 4 (/i)//i, and the rotational 
temperature, T r (fi) = A r (/i)//(/i), where I(fi) is the moment of inertia of a grain. Hereafter 
we will ignore the dependence of the normalized energy, a, on /i: this will be the case 
of particular examples discussed below. Then, as it follows from ( |30"D the translational 
temperature tends to a fixed value 

^tl^oo -> j|a (33) 

The asymptotic value of the rotational temperature is determined by the dependence of 
its moment of inertia on the mass, i.e., by the inner structure of a grain. Of interest are two 
cases. Suppose that the grain radius is independent of its mass: this will be referred to as 
a spongy grain. Then, I = \iia 2 and Eq. ([31]) results in 

T.l^oo -> -a. (34) 

Another case is a solid or, more general, a fractal grain with /x oc a D , where D is its 
fractal dimension. Since the mass density is proportional to r D ~ 3 , the moment of inertia is 

2 D 

1 = 3 ZJT^^ 2 ^' ^ 
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It should be noted that strictly speaking our derivation of kinetic equations is applicable 
to solid (D = 3) grains only. However, the rotational temperature for fractal grains turns 
out to be independent of the fractal dimension: 



3 

T r \^oo — > ~Ot. (36) 

Now we turn to evaluation of the normalized energy, a (0). Although generally the 



solution of the kinetic equation (22) is time-dependent and can deviate from the Maxwellian 
distribution, there are reasons to ignore this deviation. Suppose as an example, that the 
collision term I a in Eq. (p2|) is represented by a sum of the Boltzmann integral describing 
interatomic collisions and let's assume some source term balancing the gas loss. Then, 
with the dominating role of elastic interatomic collisions, the deviation from Maxwellian 
distribution becomes negligible. Assuming that N n is the Maxwellian distribution with the 
temperature T , we get e = 2/3T , i.e. 

T t - ^T (37) 
I T ; fractal grain 

Of interest is that the kinetic dust temperature obtained under the same approxima- 
tions but in neglecting the mass growth Q is 2 To. Therefore, the mass growth results in 
appreciable cooling of the dust component. 

At lower gas pressure the interatomic collisions are negligible and we have to take into 
account the deviation of the ambient gas distribution from Maxwellian. Suppose there is 
some bulk source of the Maxwellian gas, i.e. I a in Eq. (|22D is given by 

UP) = v*fM(p) = ^o (27rm ^ o)3/2 exp (-^) . (39) 



Assuming that the rate of the dust mass growth (|27D is smaller than the rate of gas 
creation, i.e., mno/fiUd 1, we neglect the time derivative in (p2|), that results in 



N n (p) = ^/m(p), (40) 
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where £ = J dTira 2 (M)N ( i(T). Obviously, this distribution differs from Maxwellian due to 
the accumulation of slow atoms. Evaluating the integrals we get e = |T , that is, in the 
case of high dust density 




(41) 



T 

± r 




§T ; spongy grain 
|T ; fractal grain. 



(42) 



The main results of this section may be summarized as follows. The momentum distri- 
bution of the dust component absorbing the ambient gas tends to Maxwellian. Although 
the average value of momenta are always growing, the corresponding temperatures tend to 
certain fixed values, which differ from each other. Moreover, the rotational temperature 
depends on the inner structure of dust grains. 

In ignoring the mass growth, i.e., the last term in Eq. (p5|), the dust translational tem- 
perature in the Maxwellian gas tends to 2T ||. However, due to the absence of the friction 
torque within the adopted model, there is no steady behaviour with respect to the rotational 
degrees of freedom. 

These conclusions were gathered assuming the mass dispersion is negligible. However, 
the same inferences also follow from the general non- stationary solution of Eq. fl25|) , which 
may be obtained in a similar way but is too cumbersome to be adduced here. 

III. NUMERIC SIMULATIONS 
A. An algorithm 

The simulations of the Brownian kinetics of a single grain were performed in a following 
way. The computational area was a three-dimensional cube of unit length on edge in contact 
with unbounded equilibrium gas. This contact was simulated by point atoms of unit mass, 
which were randomly injected inside the cube from all of its faces and could freely leave the 
computational area. For each atom leaving the cube, another atom with the random velocity 
was injected from the random point of the random cube face. The distribution function of 
the injected atoms was semi-Maxwellian. 
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It was verified that in the absence of dust grains the bulk distribution inside the cube was 
Maxwellian with the prescribed temperature, T . The average number of atoms depended 
on the thermal velocity; in most runs it fluctuated around 10000. Since there were no forces 
acting upon atoms, their trajectories were straight lines. 

The grain was represented by a movable sphere of a radius small compared to the cube 
edge, typically, the initial size was a = 0.01. If an atom hit the grain surface then it 
transferred a part of its momenta to the grain according to some prescribed rules, which are 
discussed below. In the case of absorbing grain, no new atom was injected into the cube 
after the collision, that resulted in some reduction of density. Since the grain was initially 
small, a very little part of atoms could actually experience the collision. The equations of 
motion for both translational and rotational degrees of freedom of the grain were solved; the 
time step was small compared to the average time between the collisions. 

It was observed that some spurious force arose when the grain approached the faces of the 
cube. To minimize the influence of this computational effect the grain was confined near the 
center of the cube with the help of the auxiliary spherically-symmetric parabolic potential 
well. The parameters of the well were chosen in such a way that even with the kinetic 
energy of about 100 T the grain could not approach the faces of the cube. Evidently, such 
a confinement should result in multiplication of the grain distribution by a Boltzmannian 
factor and could not alter the distribution over the kinetic energy. Moreover, in application 
to the dusty plasma the confinement appears in a natural way. The imposed auxiliary 
potential did not influence the motion of atoms. 

It is worth noting that although the mass of the grain could be as large as several millions 
of atom masses, its mobility was crucial for our simulations. It was the temptation to avoid 
the solution of the grain equations of motion by simply counting down the energy and 
momentum transferred to the immobile grain. However, this way led nowhere: the grain 
temperature was permanently increasing without any saturation. 

The main goal of our simulations was to accumulate enough data in order to reconstruct 
the grain distribution function over its kinetic and rotational energy. The translational or 
rotational energy axis, say, < E < 20 T , was split in a number of sub-bands (usually, 50). 
Two methods of averaging were used. First, we could trace the energy variation of a single 
grain and evaluate the time it spent in each energy sub-band. Then, these time intervals 
were summed up resulting in time-averaged distribution function. 
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Another method of averaging was the simulation of the canonical Gibbs ensemble. Ini- 
tially, the grain was situated at the center of the cube, its rotational and kinetic energies 
were chosen randomly using the random-number generators of Maxwellian distributions with 
corresponding initial kinetic, T t0 , and rotational, T r0 , temperatures. The evolution of the 
grain energies was recorded for the time period < t < t max . The obtained dependence rep- 
resented a single sample from the Gibbs ensemble. The whole procedure was repeated many 
times with varied random initial energies but the same initial temperatures. By counting 
down a number of samples in each energy sub-band for a given time instant we were able to 
reconstruct the time evolution of the energy distribution function. 

It should be noted that reconstruction of the distribution over both kinetic and rotational 
energies, f(E t , E r ) requires too many samples or too long integration time. For this reason 
we could evaluate only the distribution over kinetic, f(E t ), or rotational, f{E r ), energy 
separately. 

B. Atom-grain collisions 

As it was mentioned, for test (and fun) purposes we have used various laws of interaction 
of an atom with the grain surface. 

Specular reflections. The atom experienced the specular reflection in the grain reference 
frame: 



p = p- Ap, 

P' = P + Ap, (43) 
G' = G, 

M' = M, 

where the transferred momentum is 

/ m \ 
Ap = 2n [n ■ p n ■ P 

\ M J 

and n = (r — R)/a is a unit vector. 

Absorbing grain. The grain momenta were changed according to Eq. (0). We were able 
to use various dependencies a(M). 
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FIG. 1: Energy distributions for various kinds of atom-grain interactions. The distributions over 



kinetic energy are plotted for the cold grain (•), for the specular reflections given by Eqs. (43) 
(■) and for absorbing grain ("^). The distribution over rotational energy is shown for the case of 
absorbing grain (♦). The solid curves correspond to Maxwellian distributions with T/Tq = 2 (1), 
T/T = 1 (2) and T/T =4/5 (3). 

The "cold grain" is also described by Eqs. ([2J) but its mass remained unchanged; the same 
law of collision was used in ||. Physically, this corresponds to the diffuse scattering with 
the complete energy accomodation at the cold surface, when the net atom momentum is 
transferred to the grain and the influence of the scattered atom is negligible. 



C. Simulation results 



The simulations with specular reflections fl4~3| ) were used to check the overall operation 
of the code. Both averaging methods described above were used. It was observed that for 
the relatively heavy grain (M/m = 100) after some thousands of collisions the distribution 
function over the kinetic energy eventually tended to the Maxwellian distribution. The dust 
temperature was equal to T with the accuracy less than 1%, as expected for the system in 
the thermodynamic equilibrium. There was no relaxation in rotational degrees of freedom 
because there was no coupling between translational and rotational motions. 

The simulations with absorbing grain were performed with the solid grain, i.e. a(M) oc 
M 1 / 3 . Although the computer facilities allowed us to monitor the motion of the grain for a 
very long time, up to tens of millions collisions, it was found that no statistically significant 
result could be obtained with time- averaging method. The reason is fairly evident: with 

17 



T 

1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 

10 20 30 40 50 60 70 80 90 t 

FIG. 2: Temporal dependence of translational (solid line) and rotational (dashed line) tempera- 
tures. 

growing mass the grain motion slowed down, and it took more and more time for the grain 
to migrate from one energy sub-band to another. 

The ensemble averaging yielded more meaningful results. Several typical energy dis- 
tributions in logarithmic scale are depicted in Fig. |l|. The fit curves there correspond to 
Maxwellian distributions, i.e., F(E) oc y/Eexp(— E/T). Fig. ^| shows the temporal evo- 
lution of translational (solid line) and rotational (dashed line) temperatures. The initial 
temperatures are T t o = 0.8 To and T t q = 1.3 To. The time in Fig. |^ is measured in the units 
of characteristic time of the mass growth, r M . According to Eq. (P?P, r M = M /(47raQ J), 
where M and a are the initial grain mass and radius, respectively. The averaging was 
performed over the ensemble consisting of 1000 samples. 

There was the sharp drop in translational temperature evolution at the beginning of the 
computation. It originated due to the parabolic potential well, in which the grain was kept. 
Initially the grain was situated at the bottom of the well, so its potential energy was zero. 
Then, according to the virial theorem, half of its kinetic energy transferred to the potential 
one. Since the potential well did not influence the grain rotation, there was no rotational 
energy drop. 

Eventually the dust translational temperature tends to ~ 0.75 T , that is, near the value 
given by Eqs. (|37|). The rotational temperature is also close to the value fl3"8"|). Since there 
was only one small grain, the distribution of projectile atoms was essentially Maxwellian 
and we were unable to reproduce the solutions ( f4l|) and (|42|). 

The simulations with the "cold grain" also confirmed the analytical solution: the trans- 




18 



lational temperature tends to 2T while the rotational temperature is always growing. 

IV. CONCLUSION 

Although the simple model of the atom-grain inelastic collisions accepted in this paper 
ignores some essential processes, it demonstrates many interesting features. We confirmed 
that the translational temperature of the dust component differs from the temperature of 
the ambient gas. However, we demonstrated that the "equilibrium" temperature is highly 
sensitive to the details of the inelastic collision. In particular, taking into account the grain 
mass growth results in appreciable reduction of the dust temperature. It should be noted 
that the difference between various collision laws discussed above numerically is very small. 
The discrepancies in energy balance of the order of a fraction of a tiny mass ratio are 
accumulated and eventually result in considerable effect. It is also worth pointing out that 
the rotational temperature is sensitive even to the details of the inner structure of grain. 

The important lesson, which may be drawn from the above discussion, is that there is 
no thermodynamic equilibrium between dust and ambient gas. The statement itself is fairly 
evident since a dusty plasma is an open system. However, this indicates the inapplicability 
of the fluctuation-dissipation theorem, which is the basement of the Langevin approach to 
the theory of Brownian motion. Therefore, the problem of deducing the Langevin equation 
applicable to dusty plasmas arises. 
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